Gravitational bouncing of a quantum ball: Room for Airy's function Bi 

Haret C. Rosu 

Instituto de Fisica de la Universidad de Guanajuato, Apdo Postal E-1^3, Leon, Guanajuato, Mexico 



I apply (i) a classical version of the Ermakov-Lewis procedure and (ii) the strictly isospectral supersymmetric approach to 
the Schroedinger free fall of the bouncing ball type. In both cases, the Airy function Bi, which in general is eliminated as being 
unphysical, plays a well-defined role. Relevant plots are displayed. 
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I. INTRODUCTION 

Bouncing of a cold atomic cloud has been first ob- 
served in the laboratory in 1993 Q|. It has been shown 
that cold atoms dropped onto an "atomic mirror" can 
be used for holographic manipulation of atomic beams 
pL More recently, bouncing Bose-Einstein condensates 
have been examined in the laboratory ||. Here, we con- 
sider the toy model of a Schroedinger quantum particle 
bouncing on a perfectly reflecting surface in a linear grav- 
itational field, which is known as the quantum bouncing 
ball [QBB] problem 0). In the QBB case, one should 
solve the Schroedinger equation with the potential 



Vqbb(z) — mgz, if z > 
Vqbb(z) = oo , if z < . 



(1) 



By the scalings s = z/l g and S = E/mgl g , where 
/ , 2 \V3 

^9 = ( 2m i g ) ^ s ^ e "gravitational length" unit, the 
stationary QBB Schroedinger equation becomes dimen- 
sionless 



d 2 ifj 



S)^ = 



(2) 



The general solution is a superposition of Airy functions 
Ai(s) and Bi(s), but Airy's Bi is discarded for going to 
infinity at large s. Moreover, the perfectly reflecting 
boundary requires the wave function be zero at the ori- 
gin and therefore the physical eigenmodes are written 
as V'n(s) = N n Ai(s — S n ), where N n is the normalization 
constant and S n are the zeros of the Ai function |Q. In 
other words, a shift of the Airy's argument is performed 
placing the Airy zeros at the origin. To the best of the 
author's knowledge all the previous works in this field 
made use of only Airy function Ai of shifted argument. 
The main purpose here is to show that there are two tech- 
niques in which the Bi function could still be employed 
without leading to unphysical results. One of them is 
the Ermakov-Lewis (EL) procedure, which is presented 
in section II and the other one is the strictly isospectral 
supersymmetric (SUSY) approach enclosed in section III. 
A small conclusion section ends up the work. 



II. CLASSICAL ERMAKOV-LEWIS APPROACH 
FOR QBB 

I will use the version of the EL approach fa] that I 
introduced in previous works in collaboration fqjT Eq. (2) 
can be mapped in a known way to the canonical equations 
for a classical point particle of unit mass, generalized 
coordinate q = ip, momentum p = tp, (i.e., velocity v = 
■0), where the dot means total derivative with respects 
to s, i.e., we identify the coordinate s with the classical 
Hamiltonian time. Thus, one is led to 



dq 

q =d^ = P 



P 



S-f- 8 * 



(3) 
(4) 



These equations describe the canonical motion for a clas- 
sical point particle as derived from the time-dependent 
Hamiltonian of the inverted oscillator type 



Hd(s) = ^-(s-S)^ 



(5) 



For this classical Hamiltonian the triplet of phase-space 

2 2 

functions T\ = \, T% = pq, and T 3 = forms a dynam- 
ical Lie algebra, i.e., H c i = 2 n =i n n(s)T n (p, q), which 
is closed with respect to the Poisson bracket, namely 
{Ti,T 2 } = -2T U {T 2 ,T 3 } = -2T 3 , {T 1; T 3 } = -T 2 . Us- 
ing this algebra H c i reads 



H r 



Ti - (s - S)T 3 



(6) 



The Lewis invariant X belongs to the dynamical algebra, 
i.e., one can write T(s) = J^ r er(s)2r, and by means of 
|^ = — {I, H} one is led to the following equations for 
the functions e r (s) 



e„ = , 



(7) 



where CJj m are the structure constants of the Lie algebra 
that have been already given above. Thus, we get 



1 



ei = -2e 2 

e 2 = -(s - S)ei - e 3 
e 3 = -2(s-S)e 2 . 



(8) 



The solution of this system can be readily obtained by 
setting ei = p 2 giving e 2 — —pp and e 3 = p 2 + -L, where 
p is the solution of the Milne-Pinney (MP) equation 
0, p — (s — S)p = \. Since Pinney's note in 1950 it is 
widely known how to write p as a function of the two 
particular solutions of the corresponding parametric os- 
cillator problem. We have followed the method of Eliezer 
and Gray ||] in order to write pis) as a combination of 
Airy functions that satisfy the initial conditions as given 
by those authors. Thus, we used in all our calculations 
the following formula 



Pl (s) = Ni [(Ai(s - Sr) + Bi(s - Sr)) 2 + Bi 2 (s - Si)] 



1/2 



(9) 



i.e., we used the two Airy functions corresponding to 
the ground state. In Eq. (9), Si = (9tt/8) 2/3 and 
Ni = (8tt 2 /9) 1/6 [§. In terms of the MP solution p(s) 
the Lewis invariant reads 



2n(s) 



(PnP - p n q) 2 



2P 2 



(10) 



For example, one can check by direct calculation that 
I\ (s) = j since according to the Eliezer-Gray interpreta- 
tion the EL invariant should be ^h 2 where h is the coeffi- 
cient of the inverse cubic nonlinearity in the aforewritten 
MP equation where h = 1. 

In the EL approach the angular quantities are given 
by the following formulas 0] 



A6 d 



1 1 d . .o 



ds 



and 



A6»s = - 



' ^±{ PP )-2p 2 
ds 



ds 



(11) 



(12) 



for the dynamical and geometrical angles, respectively. 
Thus, the total angle will be 



III. STRICTLY ISOSPECTRAL SUSY 
BOUNCING BALL 



Factorizations of one-dimensional Schroedinger oper- 
ators have been first discussed in the SUSY context by 
Witten in 1981 |l0| , and are well known in the mathemat- 
ical literature in the broader sense of Darboux covariance 
of Schroedinger equations Jll| . 

In 1984, Mielnik [Q introduced a different factoriza- 
tion of the quantum harmonic oscillator based on the 
general Riccati solution here denoted by w g . As a result, 
Mielnik obtained a one-parameter family of potentials 
with exactly the same spectrum as that of the harmonic 
oscillator. Mielnik's method offers an interesting possi- 
bility to construct families of potentials strictly isospec- 
tral with respect to the initial (bosonic) one by sim- 
ply taking into account the most general superpotential 
(i.e., the general Riccati solution). Thus, in the QBB 
case one requires V+(s) = w 2 + where V + is the 
fcrmionic partner potential of Vqbb- It is easy to see 
that one particular solution to this equation is w p = w(s), 
where w(s) = —tp[/ipi is the common Witten superpoten- 
tial. One is led to consider the following Riccati equa- 
tion w 



g 



ds 



ds 



whose general solution can 
be written down as w g (s) = w p (s) + where v(s) is 
an unknown function. Using this ansatz, one obtains for 
the function v(s) the following Bernoulli equation 



dv(s) 
ds 



2v(s) w p (s) = 1, 



that has the solution 

v(s) 



Io(s) + A 
V 2 (s) 



(14) 



(15) 



The integral Io (s) = J Q S tp 2 (y) dy is a step- like function 
as one can see in Fig. 4. On the other hand, A > 
is an integration constant thereby considered as a free 
parameter, which is a measure of the contribution of the 
second linearly independent solution, i.e., the Airy Bi in 
the QBB case, as we argued elsewhere |l3| . Thus, w g (s) 
can be written as follows 



w g (s; A) = w p (s) + — ln(I (s) + A) 



d 

ds 



In 



ds 

Ms) 



Io(s) + A 



(16) 



Ati' = AH' 1 -f At) K : / ^ds . 

P 



(13) 



Plots of all these angles calculated using p\ are displayed 
in Figs 1,2,3, respectively. 



Finally, one easily gets the parametric family of poten- 
tials 



V(s; A) = w 2 (s;A) 



dw g (s; A) 
ds 



Vqbb(s 



ds 2 



ln(I (s) + A) 



2 



Vqbb(s) 



4V>i(s)V4(s) , 2^f(s) 



I (s) + A (Io(s) + A)2 



(17) 



All V(s; A) have the same SUSY partner potential V+(s) 
obtained by deleting the ground state. They may be 
considered as a sort of intermediates between the bosonic 
potential Vqbb(s) and the fermionic counterpart V+(s). 
A plot of V(s; A) is given in Fig. 5. From Eq. (16) one can 
infer the ground state wave functions for the potentials 
V(s; A) as follows 



¥>i(s;A)=N(A) 



Ms) 

Io(s) + A 



(18) 



where N(A) is a normalization factor that can be shown 
to be of the form N(A) = \/X(X + 1). The normalized 
functions ip\ and ip± are plotted in Fig. 6. 



IV. CONCLUSION 

Airy's function Bi can find a place in the physics of 
the quantum bouncing ball through two theoretical pro- 
cedures connecting the Schroedinger equation with the 
nonlinear Milne-Pinney equation and Riccati equation, 
respectively. This may help in gaining further insight 
in the problem of nonrelativistic quantum free fall. The 
Lewis angles and phases that depend on the function Bi 
through the Milne-Pinney function are important quan- 
tities provided by the Ermakov-Lewis approach that is 
here applied to a Schroedinger free fall problem for the 
first time. These quantities are similar to Berry phases 
and Hannay angles and in principle can be measured in 
quantum bouncing ball experiments. On the other hand, 
in the strictly isospectral supersymmetric method, the 
contribution of the Bi function enters through the param- 
eter A jn|. However, although the results are physically 
sound, it is still not clear what is the corresponding ex- 
perimental configuration. In other words, it is not clear 
how a strictly isospectral partner potential, such as the 
one displayed in Fig. 5 can be produced experimentally. 
For example, one may think of some particular micro- 
scopic surface effects of the atomic mirror that might be 
able to distort the interaction potential in the way the 
SUSY scheme predicts. 
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FIG. 0. Ermakov-Lewis invariant Xi(s) cf. Eq. (10) [not in 
the accepted version. 



FIG. 3. Lewis' total angle cf. Eq. (13). 
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FIG. 1. Lewis' dynamical angle cf. Eq. (11). 



FIG. 4. The integral Io (s) of the strictly isospectral SUSY 
QBB. 
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FIG. 2. Lewis' geometric angle cf. Eq. (12). 



FIG. 5. The strictly isospectral QBB gravitational potential 
for A = 1. 



FIG. 6. The normalized wave functions tpi(s) (full line) and 
ipi (s; 1) (dashed line). 
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